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In a very long Gaussian polymer on time scales shorter that the maximal relaxation time, the mean 
squared distance travelled by a tagged monomer grows as ~ t^^"^ . We analyze such sub-difTusive 
behavior in the presence of one or two absorbing boundaries and demonstrate the differences between 
this process and the sub-diffusion described by the fractional Fokker-Planck equation. In particular, 
we show that the mean absorption time of diffuser between two absorbing boundaries is finite. Our 
results restrict the form of the effective dispersion equation that may describe such sub-diffusive 
processes. 
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I. INTRODUCTION 

The stochastic fluctuations of a broad range of physi- 
cal systems [ij exhibit a behavior commonly denoted as 
anomalous diffusion. The random motion is character- 
ized by the scaling of a mean squared coordinate, which 
(when averaged over many reahzations) scales as in 
time t. For "normal" diffusion a = 1, while the cases of 
a 7^ 1 are referred to as "anomalous," with a < 1, cor- 
responding to sub-dijfusion, while a > 1 describes super- 
diffusion. The physical origin of anomalous behavior is 
usually the coupling of the particle (or some other coor- 
dinate) to many other degrees of freedom, such that its 
dynamics is the superposition of numerous other modes 
with widely distributed time scales. In principle, there is 
no reason to expect any "universality" in such anomalous 
processes, and any two situations described by the same 
exponent a may have very different characteristics. Nev- 
ertheless, certain general considerations have motivated 
approaches that encompass a large variety of cases: for a 
review see Ref. . 

Many problems related to the behavior of random 
walkers can be formulated in terms oi first passage of a 
walker, or its absorption at a boundary [3,|3,[5[. The pres- 
ence of the absorbing boundary may help to discriminate 
between different types of anomalous random walkers. 
Indeed, in the following subsections we shall demonstrate 
how the study of absorption can be used to gain better 
understanding of the complexity of anomalous behavior. 
The behavior of a normal diffuser confined by absorb- 
ing boundaries is well understood; in particular, for large 
times t the survival probability S{t) of such a diffuser 
decays exponentially. Consequently, the probability den- 
sity function (PDF) of the diffuser to be absorbed at a 
particular time Q(t) — —dS{t)/dt also exhibits an expo- 
nential decay, leading to a finite mean absorption time. 
The corresponding result in the case of sub-diffusion is 
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less clear. Only recently it was established Q (while 
building on previously known expressions 0, [zl ) that for 
one-dimensional (ID) sub-diffusion between two absorb- 
ing boundaries which is described by a particular (frac- 
tional) diffusion equation , the PDF of absorption Q{t) 
decays as a small power of t, leading to an infinite mean 
absorption time. The fractional diffusion equation used 
in this analysis [6] applies to continuous time random 
walks, which at each step have a waiting time distribu- 
tion with a long tail. 

A relatively simple and practically important case of 
sub-diffusion is the motion of a tagged monomer in a long 
polymer, whose anomalous dynamics was deduced and 
(numerically) observed by Kremer et al. Q. A polymer 
consisting of a large number N of monomers has pro- 
cesses happening on multiple length scales, ranging from 
the microscopic distance, such as separation between ad- 
jacent monomers, a, to the size of the polymer. (A mea- 
sure of the latter is the radius of gyration Rg. In a good 
solvent Rg ~ aN^ |§], with the exponent v w 0.59 in 
space dimension d = 3. The "«" sign indicates omission 
of a dimensionless prefactor of order unity. In the ab- 
sence of inter-monomer repulsion v = \/2 for any d.) To 
these length scales are associated times Tmicro ~ o? I Do, 
below which a selected monomer "does not feel" its sur- 
roundings, and Tat « i?g/DcM ~ a^N^'^'^^ / Do for how 
long it takes the polymer to diffuse its own Rg. Here 
Do denotes the diffusion constant of a single monomer, 
while the diffusion constant of the entire polymer, or its 
center of mass (CM), is Do/N . (In this discussion we dis- 
regard hydrodynamic interactions.) Very short and very 
long times correspond to normal diffusion with different 
diffusion constants. It has been shown in Refs.[8] that 
for intermediate times Tmicro < t < tn, the polymer un- 
dergoes anomalous diffusion with mean squared distance 
« a'^-^°'{Dot)°', where a = 2iy/{l -\- 2v). Note, that sub- 
diffusion occurs even in the case of the ideal polymer with 
V ^ 1/2. 

In this work we analyze the sub-diffusive motion of 
a tagged monomer which is part of an ideal (Gaussian) 
ID polymer. While the entire polymer performs diffu- 
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sive (Monte Carlo) dynamics, we record the position of 
a single monomer at the mid-point of the chain. In the 
absence of absorption, this is a simple, analytically solv- 
able problem, and the exact PDF of the monomer posi- 
tion is easily obtained. We could not extend these so- 
lutions in the presence of absorbing boundaries, and in- 
stead resorted to numerical studies. Indeed, with a single 
absorbing boundary, the PDF of the monomer position 
cannot be found using the standard method of images, 
which is the standard approach for normal diffusion and 
even some cases of sub-diffusion. We show that when 
such a particle is placed between two absorbing bound- 
aries, it has a finite mean absorption time, which scales 
(as expected) with the distance between the absorbing 
boundaries. Thus, the tagged monomer presents a sim- 
ple example of sub-diffusion whose survival probability 
differs drastically from that obtained by application of 
fractional diffusion. 

To provide the basis of comparison with anomalous dif- 
fusion, in Sec. |ll] we briefly review the behavior of a nor- 
mal diffuser in the presence of absorbing boundaries. Our 
model of the tagged monomer, and the numerical proce- 
dure used, are presented in Sec. lIIII We also present some 
numerical results confirming the expected sub-diffusive 
motion of a single monomer. In Sees . HVl and IVl we study 
the behavior of the tagged monomer in the presence of 
a single, and a pair of absorbing walls, respectively. We 
thereby demonstrate the similarities and distinctions be- 
tween our anomalous diffuser and a normal random walk. 
Notably, we stress the differences between our case and 
the solution to the fractional diffusion equation. In the 
final Sec. IVH we discuss the possible applicability of our 
results to the the translocation of a polymer through a 
membrane pore, which was in fact one of the motivations 
for this study. 



II. NORMAL DIFFUSION WITH ABSORBING 
BOUNDARIES 



The simplest model of a Brownian particle is a ran- 
dom walk (RW) on a discrete lattice, in which both the 
position of particle R and time (number of steps) t are in- 
tegers. Exact expressions for the PDF p(R, t), and many 
other properties, are readily available [lOj. A continuum 
version is the Langevin equation for the motion (diffu- 
sion) of a single particle in a solvent (in the high friction 
limit) , moving under the influence of thermal noise [Tl| : 



at a later time is 



(1) 



Here C is the friction coefficient, and the thermal 
noise satisfies {ri{t)) = (no bias) and {ri{t)r]{t')) — 
2(kBTS{t — t'). (The angular brackets, (), indicate av- 
erages over different realizations of the thermal noise.) 
Starting at i? = at t = 0, the PDF of particle position 
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where = 2Dot depends on the diffusion constant 
Do = ksT / C,. From the Langevin equation, one can also 
dire ctly construct the Fokker-Planck (diffusion) equation 
[nllil for the PDF, as 
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(Throughout this paper we consider one-dimensional mo- 
tion; the generalization to higher dimensions is straight- 
forward.) 

Consider a diffusing particle starting at the origin and 
reaching position R at time t without ever touching an ab- 
sorbing boundary at i?a > 0. The solution of both contin- 
uous and discrete versions of this problem have been de- 
scribed in detail by Chandrasekhar [13]. It can be shown 
that the PDF for the random walker must vanish on the 
absorbing boundary. Since the diffusion equation is lin- 
ear, this boundary condition can be satisfied by superpos- 
ing the PDF of a free particle (Eq. ([2])), and one starting 
at a reflected image as P(i?, t) = P(i?, t) - P{R-2R^, t). 
The survival probability S{t) = ^^^P{R,t) reduced by 
absorption [4], and for large times the absorption PDF 
decays as Q{t) oc t~^/'^. This PDF has diverging mean, 
since the particle can drift infinitely far in the direction 
opposite the wall. It should be noted that the image 
method is specifically suited to random walkers perform- 
ing independent unit steps] it fails for the sub-diffusive 
walkers considered in this paper, and also for long-range 
hops of super-diffusive motion [lj| . 

If the ID diffusing particle is confined by two absorb- 
ing boundaries, one can still use superposition by the 
method of images to create a solution. However, in order 
to satisfy both boundary conditions, an infinite set of im- 
ages is necessary. A more convenient answer is obtained 
by expanding the solution in terms of the eigenfunctions 
of the diffusion equation ([3]). For a particle enclosed by 
absorbing boundaries at i?ai = and _Ra2 = this gives 
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where {An} depend on the initial conditions. Note that 
at long times P{R,t) ~ Ai sin(7rP/L)e-(''/^)'^°*, i.e. 
the PDF has a simple sinusoidal shape with zeroes on 
the boundaries. The survival probability is S{t) oc e^*/'^, 
where the characteristic decay time r = /tt'^Do is of 
the order of time the particle needs to diffuse over the 
length of the interval. The PDF for absorption, Q{t), 
also decays with the same time constant. 

Anomalous diffusion can in principle have a myriad of 
distinct causes. An extensively studied case corresponds 
to the so-called continuous time random walks for which 



the waiting time between successive steps is taken from 
a broad distribution, with power law tails and a diverg- 
ing mean. The interest in such processes originated in 
studies of diffusion in semi-conductors [TB'l, but they even- 
tually became a prototype of anomalous diffusion. The 
fractional diffusion equation (which involves an integral 
operator) was developed to describe the evolution of the 
PDF for such walkers [la . KLTII . and explicit solutions are 
now available 0, [13, oEMi- Unlike Eq. © the solu- 
tion to these equations is not smooth, but has a cusp at 
the origin. Another interesting feature is the behavior of 
the absorption PDF for anomalous diffusers between two 
absorbing boundaries: it has been shown Q, by careful 
analysis of the solutions [3, that for large times, Q{t) 
decays as t~^^~^"\ leading, for a < 1, to a diverging mean 
absorption time! 



III. THE MODEL AND NUMERICAL 
PROCEDURE 

A monomer in a polymer undergoes anomalous dif- 
fusion even in the trivial case of a phantom chain with 
no interactions, for which v = 1/2 in any d. For this 
value of v, monomer fluctuations are governed by an ex- 
ponent a = 1/2, and thus exhibit sub-diffusion. Since 
many properties of phantom polymers can be calculated 
exactly, this presents an excellent model for the study of 
sub-diffusion. For simplicity, we shall restrict ourselves 
to the one-dimensional situation; generalization to higher 
dimensions is straightforward. 

After coarse-graining, a sufficiently long flexible poly- 
mer can be represented by effective monomers connected 
to their nearest neighbors by harmonic potentials (Gaus- 
sian springs) 0, [lOl . Thus, the Hamiltonian for a chain 
of A^-monomers is 

K 

H^-Y^iRn+l-Rn)'. (5) 
^ n=l 

The distribution of the distance between two adjacent 
(along the chain) monomers at a temperature T is gov- 
erned by the Boltzmann factor exp[— /3if(i?„+i— i?„)^/2], 
with P — \/{kBT) and ks is the Boltzmann con- 
stant. The mean squared separation between adjacent 
monomers is a? = ksT/K, while the mean squared 
radius of gyration is i?^ = i7v(l - l/N'^)kBT/K « 
\NkBT/K. 

Theoretical treatment of the polymer described above 
requires solution of N coupled Langevin equations. How- 
ever, the problem becomes particularly simple if we de- 
scribe the configurations using Rouse modes [2(| 

1 ^ 1 

where q = pn/N, and p = 0, 1, . . . A^ — 1. In terms of 
Uq, Langevin equations decouple, and each Rouse mode 
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FIG. 1: (Color online) PDF for the position of the central 
monomer in a Gaussian chain of 129 monomers at times t = 
2000, 4000, ■ ■ ■ , 20000 (from narrowest to the widest graph). 
The distributions were obtained from 100,000 runs using bin 
size AR = 0.1. The continuous lines represent normalized 
Gaussian fits to the distributions. 



can be viewed as an independent "particle" moving in a 
harmonic potential whose strength depends on q. The 
PDF of every Uq is Gaussian. Conversely, the position 
of each monomer can be viewed a linear combination of 
UqS (inverse of Eq. ©). Since the linear combination of 
Gaussian variables is a Gaussian variable, we are assured 
that each monomer is described exactly by a Gaussian 
PDF, and the theoretical study is reduced to evaluation 
of the mean and variance of that distribution (see later). 
In the presence of absorption, such treatment is not pos- 
sible. 

We measured distances in dimensionless units, i.e. 
multiplied by -y/ K/ksT. In these units the the root 
mean square separation of adjacent monomers is a = 1, 
and Rg « A^/6. We used diffusive (Rouse) dynamics to 
evolve the system in time, i.e. the monomers were moved 
using standard Monte Carlo (MC) moves. An elementary 
MC move consists of randomly picking one monomer and 
attempting to increment its position by SR chosen uni- 
formly from the interval (-1,-1-1), in dimensionless units. 
The change in the Boltzmann weight factor controls the 
probabilistic decision of whether the move is accepted. A^ 
elementary move attempts are defined as one MC time 
unit. The mean squared displacement of a monomer in 
a single move determines the diffusion constant Do; with 
the above choice of step size we had Do = 0.10. 

For simulations we chose polymers of odd lengths A^ — 
2^ + 1, with £ = 1, 2, • • • , 10, i.e. A^ = 3, 5, 9, • • • , 1025. 
While all the monomers moved during the simulation, we 
followed only the position of the central monomer num- 
bered c = 2^~^ -I- 1. For each case, 100,000 independent 
simulations were performed to ensure reliable averages. 

As an example. Fig. [1] depicts the PDF of the position 
of a central monomer (c — 65) in a Gaussian polymer of 
A^ — 129 monomers. At t = all monomers were located 
at the coordinate origin. As the configuration of the poly- 
mer evolved in time, the position of the 65th monomer 



10^ 



-20 -10 

position 



FIG. 2: (Color online) Logarithmic plot of the variance of the 
probability distributions depicted at Fig. [T] as a function of 
time. 



was recorded. Repeating the process 100,000 tiraes pro- 
duced the distributions shown if Fig. [TJ Note the excel- 
lent (single parameter) fit of the normalized Gaussian to 
the actual graphs. Indeed, as explained in the discus- 
sion following Eq. ([6]) the PDF of the monomer must be 
Gaussian at all times. It should be noted, that this shape 
is the same as in the case of the normal diffusion, and 
significantly differs from solutions of sub-diffusive frac- 
tional diffusion equations which contain a cusp at the 
origin (see, Ref. L2.])- 

While the shapes of the graphs in Fig. [1] are not 
anomalous, the time-dependence of their variance is. For 
times shorter than the longest relaxation time r^r = 
Rg/2DcM ~ a'^N^/12Do, which in the above case be- 
comes Ti29 — 1.4 X lO'*, the variance of the distribution 
grows as t^/^, while for times longer than rjv it is linear 
in t. This result can be demonstrated analytically, since 
the variance of the particle position can be expressed as 
a sum of variances of Rouse modes. (Analogous calcu- 
lation for a fluctuating line (or surface) can be found in 
Ref. [m.) Figure [2] depicts the dependence of the vari- 
ance on t. While all the points are within a half-decade 
from the crossover point, one can clearly discern the two 
types of behavior: the slope of the straight line through 
the first 4 points is 0.52, very close to the expected 1/2, 
and gradually increases to the right of the graph. 



IV. A SINGLE ABSORBING BOUNDARY 

Let us now introduce absorption into the problem. We 
assume that at i = all monomers are located at i? = 0, 
and an absorbing boundary is placed at i?a = 8, i.e. when 
the central particle reaches this point it is absorbed and 
the diffusion process ends. It should be stressed that 
other monomers of the polymer do not feel the absorbing 
boundary; their sole function is to generate anomalous 
diffusion of the tagged particle. 

We begin with a very short polymer with TV = 3, whose 
radius of gyration is significantly shorter than i?a- More 



FIG. 3: (Color online) PDF of the position of the central 
particle in a polymer with = 3 monomers with an absorbing 
wall at _Ra = 8. The graphs correspond (narrow to broad) 
to t = 200, 400, ■ ■ • 1000, and solid lines represent fits to the 
difference of two identical Gaussians centered at i? = and 
R — 16. The distributions are obtained from 100,000 runs, 
and the bin size is AR — 0.1. 



importantly, its maximal relaxation time T3 w 8 is sig- 
nificantly shorter than the time (about 10"^) for the CM 
of the polymer to diffuse the distance from the origin to 
the absorbing boundary. Therefore, at the time-scales 
at which the particle can be absorbed, the motion of the 
tagged monomer is indistinguishable from that of the CM 
of the polymer. Consequently, the problem of absorption 
of the central monomer should be indistinguishable from 
that of a normal diffuser with diffusion constant Do/3. 

Figure [3] depicts the observed PDF of the position of 
the central monomer of a polymer with A'^ = 3 at var- 
ious times. The area under the graphs decreases with 
time due to absorption. The shapes are the same as ex- 
pected for normal single particle diffusion: the solid lines 
are (single parameter) fits to a difference between two 
Gaussians centered at image points i? = and R = 2i?a. 
The excellent fits demonstrate that normal diffusion well 
describes the absorption for such small N. Moreover, 
the variances of the Gaussians fits increase linearly with 
time, with a prefactor corresponding to cr^(i) = 2DcMt, 
in which Dqm was calculated independently. 

This behavior changes radically when N becomes 
large. Already for N = 129 all resemblance to reg- 
ular diffusion vanishes. The maximal relaxation time 
^129 = 1.4 X 10^ is of the same order as the time re- 
quired for the CM of the polymer to diffuse the distance 
to the absorbing boundary (about 4 x 10^), and Rg of 
the polymer is of the order of the distance to the absorb- 
ing boundary. The PDF depicted in Fig. 2] cannot be 
fitted by the difference between two Gaussians at image 
points. In particular, the PDF is not linear close to the 
boundary, but appears to vanish quadratically. Thus the 
qualitative behavior changes drastically on going from 
regular diffusion for small A^ to anomalous diffusion at 
large A^. 
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FIG. 4: (Color online) PDF of tlie position of central particle 
in a polymer with A*' = 129 monomers with an absorbing 
wall at i?a ~ 8. The graphs correspond (narrow to broad) 
to times t = 2000, 4000, ■ ■ ■ , 20000. The distributions were 
obtained from 100,000 runs using bin size A_R — 0.1. 



V. TWO ABSORBING BOUNDARIES 



FIG. 5: PDF of the absorption time (in Monte Carlo time 
units) of the central monomer in a Gaussian polymer of length 
A^. The tagged monomer can be absorbed at one of two 
boundaries at i?a = ±8. The plots (from left to right, solid 
lines) correspond to N=3, 9, 33, 129, and A'^ = 513 (dots), 
and the histogram was calculated from 100,000 independent 
runs, with bin size At = 200. The graphs for N = 257, 1025 
(not shown), and for A*' = 129, are virtually indistinguishable. 



We next consider absorption of the central monomer 
by two boundaries located at i?a = ±8, for N rang- 
ing from 3 to 1025. Figure depicts on a semiloga- 
rithmic scale the PDF of the absorption Q{t) for sev- 
eral N. For small polymers the curves are indistinguish- 
able from that of a normal random walker with diffusion 
constant Dqm = Dq/N, which can be calculated from 
Eq. ([5]), with {An} selected to correspond to the ini- 
tial state (P(i?,0) = S{R)). After integrating P{R,t) 
over R to obtain S{t), we get Q{t) = —dS/dt. As ex- 
plained in Sec. |lTl for large times Q{t) decays exponen- 
tially with a time constant r « R'^/Dcm = NR^/Dq. 
The mean absorption time is of the same order of mag- 
nitude. The dependence of the mean time, and of the 
time constant for decay, is depicted by squares and cir- 
cles, respectively in Fig. [51 Note that when N becomes 
large enough, so that the longest relaxation time of the 
polymer exceeds the typical time it takes for a particle 
to travel the distance between the absorbing boundaries, 
Q{t) becomes independent of N . Indeed all the graphs 
for TV = 129, 257, 513, 1025 coincide with each other. The 
long time behavior remains an exponential decay, as can 
be seen from the straight lines on the semi-logarithmic 
plot. These curves thus depict true anomalous diffusion 
in the "infinite- iV limit," and the corresponding expo- 
nential decay time constant scales as the time it takes to 
cover the interval by sub-diffusion, i.e. r « R'^/{a^Do). 

So far, we reported on simulations in which at time 
t = 0, the entire polymer is located at the origin, i.e. 
Rn = for all n. This is a particularly convenient choice 
for analytical calculations, since all Rouse modes vanish 
at t — and their mean values (averaged over realiza- 
tions of the noise) remain zero at all times. In any case, 
we know that the initial value of each Rouse mode will 
be forgotten after one relaxation time of that mode. One 
may consider a different case, where at t = the polymer 




FIG. 6; (Color online) The mean first passage (absorption) 
time (square), and the decay time-constant of the absorp- 
tion PDF (circle), as a function of polymer size N, when all 
monomers were located at the origin at t = 0. Diamonds and 
Xs show the same quantities when at time t — only the cen- 
tral monomer is at the origin, while the remaining monomers 
are in a typical equilibrium position. The data are obtained 
from 100,000 independent runs. The statistical error bars in 
the mean times (approximately 0.3%), and the estimated sys- 
tematic errors in the decay time-constants (less than 5%) are 
much smaller than symbol sizes. 



assumes a randomly selected equilibrium configuration. 
Thus, in addition to averaging our results over different 
realizations, we also need to average over the starting 
configurations. This initial condition appears more nat- 
ural since the time t = is not special. In any case, 
we find that the differences between the two procedures 
are rather small. For small N we cannot expect much 
difference, because by the time the polymer reaches the 
absorbing boundary it is equilibrated in any case. The 
results for mean absorption time and decay time con- 
stant are depicted by diamonds and Xs, respectively, in 
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FIG. 7: (Color online) PDF of the central monomer of a 129- 
monomer Gaussian polymer absorbed at boundaries a dis- 
tance 8 from the initial position on both sides. The curve 
with the larger mean corresponds to the case when at t = 
all Rn ~ 0, while the other plot corresponds to starting config- 
uration where Rc = while the remaining monomers are equi- 
librated. This plot (histogram) was calculated from 100,000 
independent samples; bin size At = 200. 



Fig. [HI We see that the new procedure gives slightly 
shorter mean absorption times and essentially the same 
decay times. Figure [7] depicts the PDF of absorption 
times for N = 128. It seems that at short times the ran- 
dom starting point diffuser moves slightly faster leading 
to shorter mean times, but for large times both cases are 
characterized by the same decay constant. 

With its exponential decay the long time behavior of 
a particle between absorbing boundaries more resembles 
normal diffusion, although the time scales have to be de- 
termined using anomalous diffusion arguments. Never- 
theless, the PDF of the unabsorbed monomer at long 
times does not resemble that of a normal diffuser. We 
studied the PDF of the positions of surviving particles in 
100,000 independent runs for a polymer with N = 129. 
Naturally, as the time increases the probability of not be- 
ing absorbed decreases. (The decrease in probability also 
means that for large t the PDF was derived from samples 
significantly smaller than 100,000, and consequently the 
statistical accuracy of the results decreased.) To enable 
a convenient comparison between the PDFs at various 
times we normalized them to 1. In the results depicted 
in Fig. m the PDFs of the particle position were recorded 
at different times, all of the order of mean absorption 
time. Superficially these results resemble regular diffu- 
sion. In analogy to Eq. (|4]), it appears as if at very long 
times only a slowest "eigenmode" survives and the PDF 
decays as ^'(i?/i?a)e~*^^, with an eigenvalue related to r. 
The eigenfunction depicted in Fig. [51 appears to be 
universal, although specific to our form of sub-diffusion, 
while r scales as i?f and is independent of N. 

In the case of a regular diffusion (in the scaled variable 
X = R/Ra) we have '^{x) oc cos{nx/2), i.e. the func- 
tion vanishes at the boundaries {x = ±1) with a finite 
slope. By contrast, the results depicted in Fig. [5] suggest 



FIG. 8: (Color online) Probability density of the position 
of the central monomer of a 129-monomer polymer, between 
absorbing walls located at Ra = ±8. The curves are ob- 
tained at times t = 2000,4000,6000,8000,10000, and are 
the result of 100,000 independent simulations collected in 
bins of size AR = 0.1. To remove the effect of absorption 
all curves have been normalized, such that their integral is 
unity. The t — 2000 curve is relatively smooth, while that of 
t = 10000 is very noisy because only in a small portion of runs 
the tagged particle was not absorbed, and consequently the 
sample size is small. Smooth lines represent two suggested 
eigenfunctions '^{xR^) / Ra,: the dashed line corresponds to 
\l/(a;) = (1 -I- cos(7rx))/2 (see text), while the solid line rep- 
resents the normalized eigenfunction given by Eq. (A4) in 
Ref. H. 

a vanishing slope at the boundary. In fact, an attempt 
to fit the function by a few terms of the Fourier series 
6o + ^1 cos(7ra;) -f . . . gives &o ~ while the coefficients of 
higher Fourier cornponents are by an order of magnitude 
smaller. In Ref. [1J| the fractional Laplacian operator 
was examined in a bounded domain. Their particular 
implementation of boundary conditions enabled calcula- 
tion of the eigenfunction ^'(x) for various values of the 
fractional order. In particular, the explicit expression for 
the case corresponding to sub-diffusion with a = 1/2 in 
our notation is given in Eq. (A4) of Ref. fli*! with a = 4 
in their notation. While this (normalized) function, de- 
picted by a smooth solid line in Fig. [8l qualitatively re- 
sembles the numerical curves, it does not provide a quan- 
titative fit. This makes the fractional Laplacian operator 
a somewhat unlikely candidate for describing the long- 
time behavior of our diffuser. 



VI. DISCUSSION 

In this work we concentrated on an extremely simple, 
and yet non-trivial model of sub-diffusion. The Gaus- 
sian nature allows analytic calculation of some proper- 
ties, such as the probability distribution of freely moving 
particles; but the absorption properties were studied nu- 
merically. We believe that similar results should apply 
to self-avoiding polymers, although separation into inde- 
pendent Rouse modes is no longer possible, and probably 
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not much can be done beyond simple scaling arguments. 
Within our simple model, we find that the process of 
absorption is quite different from other simplified sub- 
diffusive processes in the literature. 

Anomalous diffusion of a monomer has several fea- 
tures resembling the translocation of a long polymer 
through a narrow pore in a membrane. This process 
has been extensively studied experimentally during the 
last decade [H, [23, In the theoretical descrip- 

tion of translocation, a single variable s rep resenting the 
monomer number at the pore [2^, [2^, l27l [28| indicates 
how much of the polymer has passed to the other side. 
If the translocation process is very slow, the mean force 
acting on the monomer in the hole can be determined 
from a simple calculation of entropy, and the translo- 
cation problem is reduced to the escape of a 'particle' 
(the translocation coordinate) over a potential barrier. 
Such theories produce qualitative understanding of ex- 
perimental results 29]. However, if the process is not 
slow enough, compared to the relaxation times of Rouse 
modes, then its dynamics is more complicated. Succes- 
sive steps of the reaction coordinate are then correlated 
in a manner closely resembling the correlations between 
steps of a tagged monomer in a polymer. In Ref. [S^ , it 
was numerically verified that s indeed undergoes anoma- 
lous diffusion in the ID "space of monomer numbers." 
It was further argued that the relaxation of the polymer 
constrains the translocation process and consequently de- 
termines the translocation time. Such behavior closely 
relates the translocation process to the anomalous diffu- 
sion of a single monomer. 

In the last few years significant progress has been made 
in the theoretical modelling of the translocation process. 
On one hand, short time behavior has been modelled in 
great detail |3l|, and on the other hand scaling consid- 
eration of the long time behavior have been extended 
to include hydrodynamic interactions fs^, [s^ . Recently 
Grosberg et al. [34| developed an intuitive scaling pic- 
ture of polymer translocation under the influence of a 
force. (See also Ref. [35!|.) A variety of scaling regimes 



with force applied to the end-point or at the pore have 
been investigated numerically in some detail [36|. Some 
recent studies [13, HI] suggest that the translocation pro- 
cess maybe even slower than dictated by the relaxation 
of the Rouse modes. If so, this would weaken the analogy 
between the translocation and the anomalous diffusion of 
a monomer. (The accuracy of these claims is questioned 
in further work [39].) 

To the extent that one may draw an analogy between 
translocation and anomalous diffusion in the presence of 
absorbing boundaries, one may inquire whether the mean 
translocation time is finite. Reference argues that 
translocation may be described by a fractional diffusion 
equation, and consequently require an infinite mean time, 
as found in the solutions of such equation Q. A similar 
point is made in Ref. ,38] , where a detailed study of the 
PDF of translocation times is fitted to a slowly decaying 
function for large times. However, direct (experimen- 
tal and numerical) measurements appear to indicate well 
defined average translocation times. Our results offer a 
model where absorption times of an anomalous diffuser 
arc finite. Clearly more detailed studies of translocation 
are needed to resolve this question. 

In this work we studied in detail the sub-diffusion of a 
tagged monomer in a Gaussian polymer. In the absence 
of absorption all properties can be derived analytically. 
However, upon inclusion of absorbing walls, we had to re- 
sort to numerical simulations. While we can characterize 
all properties of the numerical results, we are still miss- 
ing an equation that can describe the evolution of the 
PDF of the position of the tagged particle. In fact, the 
numerical results exclude several simple forms for such 
an equation. 
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